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Abstract Some features of the Bianchi type-I uni- 
verses in the presence of a fluid that wields an anisotropic! 

equation of state (EoS) parameter are discussed in the 
context of general relativity. The models that exhibit 
de Sitter volumetric expansion due to the constant ef- 
fective energy density (the sum of the energy density 
of the fluid and the anisotropy energy density) are of 
particular interest. We also introduce two locally ro- 
tationally symmetric models, which exhibit de Sitter 
volumetric expansion in the presence of a hypothetical 
fluid that has been obtained by minimally altering the 
conventional vacuum energy. In the first model, the 
directional EoS parameter on the x axis is assumed to 
be -1, while the ones on the other axes and the energy 
density of the fluid are allowed to be functions of time. 
In the second model, the energy density of the fluid 
is assumed to be constant, while the directional EoS 
parameters are allowed to be functions of time. 

Keywords Bianchi type-I; de Sitter expansion; cosmo- 
logical constant; anisotropic fluid 



1 Introduction 

Inhomogeneous and anisotropic universes have particu- 
larly been of interest in mathematical cosmology, rather 
than observational cosmology. This situation can be 
related to two main ingredients of modern cosmology: 
(i) The present universe had very well been described 
with Friedmann-Lemaitre models that are based on the 
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spatially homogeneous and isotropic Robertson- Walker 
space-time, (ii) According to the inflationary paradigm 
the universe should have achieved an almost isotropic 
and homogeneo us geometry at the end of the infla- 
tionary era (see iLindei 2008 for a review of inflation- 
ary cosmology). However, for a realistic cosmologi- 
cal model one should consider spatially inhomogeneous 
and anisotropic space-times and then show whether 
they can evolve to the observed amount of homogene- 
ity and isotropy. Because of the analytical difficulties 
in studying the inhomogeneous models, as a first step 
toward a more realistic model one may consider Bianchi 
type cosmological models, which form a large and al- 
most complete class of relativistic cosmological models, 
which are homogeneous but not necessarily isotropic. 
The only spatially homogeneous but anisotropic mod- 
els other than Bianchi type models ar e the Kantowski- 
Sachs locally symmetric family. See Ellis fc van Elst 
or generali zed, particu larly anisotropic, cosmo- 
logical models and Ellis! ( 2006 ) for a concise review on 
Bianchi type models. 

In principle, once the metric is generalized to Bianchi 
types, the equation of state (EoS) parameter of the fluid 
can also be generalized in a way conveniently to wield 
anisotropy with the considered metric. In such models, 
where both the metric and EoS parameter of the fluid 
are allowed to exhibit an anisotropic character, the uni- 
verse can exhibit non-trivial isotropization histories and 
it can be examined whether the metric and/or the EoS 
parameter of the fluid evolve toward isotropy. Thus, 
one gets the opportunity of constructing even more re- 
alistic models than the Bianchi type models where the 
fluid that fills the universe is assumed to be isotropic 
from the beginning. 

Moreover, in recent years Bianchi universes have 
been gaining an increasing interest o f observational 
cosmology, since the WMAP data (jHinshaw et al. 
20031 l2007i l2009h seem to require an addition to the 
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standard cosmological model with positive cosmolog- 
ical constant that resembles t he Bianchi morphology 
(seeljaffe et al 1120051 l2006allbl : ICam panelli et al.l l2006, 
l2007t iHoftuft et al.ll2009t) . According to this, the uni- 
verse should have achieved a slightly anisotropic spa- 
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1994) and that might be indicating a non-trivial 
isotropization history of the universe due to the pres- 
ence of an anisotropic energy source. We may talk 
about two main classes of such models; according to 
whether this anisotropization occurs at an early time 
or at late times of the universe. The former class can 
be related with the inflaton field, which drives the infla- 
tion, while the latter one can be related with the dark 
energy (DE) , which drives the late time acceleration of 
the universe. 

In the context of the former class, the generic in- 
flationary models can be modified in a way to end 
inflation with a slightly an is otrop i c spat ial geometry; 
see, e.g., ICampanelli et all f|2006l l2007t) . A key in- 
gredient of the inflationary models is an inflaton field, 
which generates an expansion like the one generated 
by the cosmological constant. In these models, space, 
and thus the inflaton field, is assumed to be homo- 
geneous and isotropic from the beginning. The pos- 
sibility of inflation in anisotropic space - time has first 
been investigated by iBarrow fc Turnerl (| 1 9 8 lh . Some 
authors have considered Bianchi t ype space-times in 
the pr e sence of scalar fields , e.g., iFeinstein fc Ibanez 
( 19931 ); lAguirreeabiria et al.l jl993l ). iFordl (|l989h dis 
cussed the possibility of inflationary models in which 
inflation is driven by a vector field, which gives rise 
to an anisotropic EoS parameter, rather than a scalar 
field, for the first time. However, it was suffering from 
the fin e-tuning problem. Recently, iKoivisto fc Mota 
( 2008d l have considered several new classes of viable 
vector field as alternatives to the inflaton within in the 



Bianchi type-I framework. iGolovnev et al.l (|2008l ) have 
constructed a successful model, which either could give 
a completely isotropic universe or slightly anisotropic 
universe at the end of the inflation. 

In the context of the latter class, the isotropy of 
space that has been achieved in the inflationary era 
can be distorted during the contemporary accelera- 
tion of the universe by modifying th e DE in a way to 
wield a nisot ropic pressure, se e, e.g., IKoivisto fc Mota 
( 2008allb1 1cl); iRodriguesl (|2008l ). It is well known that 
the inclusion of a positive cosmological constant, which 
is mathematically equivalent to the vacuum energy with 
p = —p and which is the simplest candidate for the DE, 
may lead to the explanation of this observed accelera- 
tion of the universe. 



Note that the vacuum energy (the cosmological con- 
stant) is coming into question in both early and late 
time acceleration of the universe. The solution of Ein- 
stein's field equations in the presence of a positive cos- 
mological constant for a homogeneous and isotropic 
model results in an exponentially expa nding universe, 
which is known as the de Sitter model (jGron fc Hervik 
2007). Thus, de Sitter and de Sitter like expansions, 



hence the cosmological constant and its alterations, are 
of particular in terest for contemporary cosmology. 
Waldl (<1983h investigated the asymptotic behavior 



of initially expanding homogeneous cosmological mod- 
els with a positive cosmological constant and showed 
that such models, of all Bianchi types (only type- 
IX under some conditions), ex ponent i ally e volve to- 
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ward the de Sitter solution, 
gated the expansion anisotropy during the inflationary 
era for a vacuum Bianchi type-I universe with a non- 
vanishing positive cosmological constant, and discussed 
the aniso t ropic generalization of the de Sitter solution. 



Beeshaml (|1994h showed that Bianchi type-I models in 
the presence of a positive cosmological constant cannot 
be of the pure de Sitter type, u nless the gravitationa l 
"constant" takes negative values. iKalligas et al.l (|1994l ) 
showed that de Sitter inflation is allowed at least at late 
times of the univers e within the Bianchi type-I frame- 
work. lArbabl (| 1997T) constructed a Bianchi type-I cos- 
mological model, in the presence of variable G and A 
and bulk viscos i ty, tha t exhibits de Sitter expansion. 



Kumar fc Singhl (|2007l) gave an exponentially expand- 
ing Bianchi type-I model in the presence o f a hypo- 
thetical perfect fluid. Akarsu fc Kilind ( 2010f ) presented 
a locally rotationally symmetric (LRS) Bianchi type-I 
model that exhibits de Sitter volumetric expansion in 
a mixture of a perfect fluid and a fluid that wields a 
special form of a dynamical and anisotropic EoS pa- 
rameter. 



B arrow! (|1997f ) presented an analysis of the cosmo- 
logical evolution of matter sources that possess small 
anisotropic pressures (electric and magnetic fields, colli- 
sionless relativistic particles, gravitons, anti-symmetric 
axion fields in low-energy string cosmologies, spatial 
curvature anisotropics and stresses arising from sim- 
ple topological defects), and he discussed the effects of 
inflation on the anisotropy of the pressures. 

Considering the above discussion, one may first think 
to distort the isotropy of the vacuum energy so as to 
obtain anisotropy in the geometry of the universe in 
general relativity. Thus, in this study, Bianchi type-I 
models that exhibit de Sitter volumetric expansion in 
the presence of a vacuum energy that has minimally 
been altered in a way to wield anisotropic EoS param- 
eter will be of our particular interest. We first, in Sect. 
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2, discuss some features of the Bianchi type-I universe 
filled with a fluid that wields an anisotropic EoS param- 
eter. We define the effective energy density as the sum 
of the energy density of the fluid and the anisotropy en- 
ergy density, and specialize the discussion to the models 
that exhibit de Sitter volumetric expansion due to the 
constant effective energy density. We then, in Sect. 3, 
examine the solutions for a LRS Bianchi type-I space- 
time, where the universe exhibits de Sitter volumetric 
expansion in the presence of a hypothetical fluid, that 
has been obtained by minimally altering the conven- 
tional vacuum energy, and we present two exact mod- 
els. In the first model, in Sect. 3.1, the directional EoS 
parameter on the x axis is assumed to be -1, while the 
ones on the other axes and the energy density of the 
fluid are allowed to be functions of time. In the sec- 
ond model, in Sect. 3.2, the energy density of the fluid 
is assumed to be constant, while the directional EoS 
parameters are allowed to be functions of time. 



2 A general discussion on the dynamics of the 
models in the presence of an anisotropic fluid 

2.1 The field equations in the presence of anisotropic 
fluid 

We consider the diagonal metric of the spatially flat, ho- 
mogeneous but anisotropic Bianchi type-I space-time, 



ds 2 



- dt 2 - A{tfdx 2 - B(t) 2 dy 2 - C(t) 2 dz 2 



(1) 



where A(t), B(t) and C(t) are the directional scale fac- 
tors, functions of the cosmic time, t. 

Within the framework of the metric given by (1), 
the energy-momentum tensor of a fluid can be written, 
most generally, in anisotropic diagonal form as follows: 



T„" =diag[To ,T 1 1 ,T 2 2 ,r 3 3 ]. 



(2) 



The Einstein field equations, in natural units (8ttG = 
1 and c = 1), are 



(3) 



where g^ v is the metric tensor; g^u^u" = 1 = 
(1,0,0,0) is the four-velocity vector); i? M „ is the Ricci 
tensor; R is the Ricci scalar, T^ u is the energy- 
momentum tensor. 

In a comoving coordinate system, the Einstein field 
equations (3), for the Bianchi type-I space-time (1) in 
case of the anisotropic energy-momentum tensor given 
by (2), read 
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where an overdot denotes d/dt. 

The directional Hubble parameters, which determine 
the expansion rates of the universe in the directions of 
x, y and z, for the metric given in (1) can be defined as 



H x — ~Ti 

A' 



H„ = — and 
y B 



H z =-, 



(8) 



respectively, and then the mean Hubble parameter, 
which determines the volumetric expansion rate of the 
universe, can be given as 



H 



1 V 
3V 



H T + H,, 



H 7 



(9) 



where V = ABC is the volume of the universe. 

2.2 The anisotropy of the expansion 

The anisotropy of the expansion can be parametrized 
by using the directional Hubble parameters (8) and the 
mean Hubble parameter (9) of the expansion, 



-At 



Hj — H 
H 



(10) 



where Hi (i=l,2,3) represents the directional Hubble 
parameters in the directions of x, y and z, respec- 
tively. A = corresponds to isotropic expansion. The 
Bianchi type-I universe approaches spatial isotropy if 
A -> 0, V -> +00 and T 00 > as t -t +oo (see 
Collins fc Hawking] 1973 for details). 

The anisotropy parameter of the expansion can also 
be given in terms of the components of the energy- 
momentum tensor (2) and the mean Hubble parameter 
(9) by using the evolution equations (5)-(7); 



9H 2 ^ 



1 2 



A,+ {TC-T^Vdt 



V 2 and i > j,\ 



(11) 



where the Xj are real constants. One can immediately 
see that the above equation is reduced to the one given 



bv lGronl(|l985h : 



_ Ai + A 2 + A 3 2 
for an isotropic fluid, i.e., where Ti 1 = T 2 2 = T 3 3 



(12) 



2.3 Anisotropy energy density in the presence of an 
anisotropic fluid 

The energy density of the fluid p, which corresponds to 
To , can immediately be written in terms of the mean 
Hubble parameter H and the anisotropy parameter of 
the expansion A by using the constraint equation (4), 



T U = 3H 2 1 



(13) 



This is equivalent to the generalized Friedmann equa- 
tion f or a Bianchi type- 1 space-time (see Ellis fc van Elstfr 
1999 and iBarrowl 11995 for the generalized Friedmann 
equation). According to this, for a given value of the 
mean Hubble parameter H in Bianchi type-I space- 
time, the anisotropy of the expansion lowers down the 
energy density; the highest energy density is achieved 
in case of isotropic expansion (i.e., A = 0) and the 
anisotropy of the expansion cannot be arbitrary be- 
cause the condition to observe p > for a comoving 
observer is A < 2. 

By moving from (13), the anisotropy energy density 
(namely, the energy density that associated with the 
anisotropy of the expansion) may be defined as follows: 



Pfi EE ^A. 

Then, using (11) in this definition, we obtain 



(14) 



1 3 

- y 



1 2 



A, + / (Tt* - T/)Vdt 



V~ 



and i > j, 
(15) 



for the anisotropy energy density. One can observe that 
the anisotropy of the fluid also contributes to the de- 
termination of the p/3 via the integral term. In case of 
an isotropic fluid the integral term vanishes and pp is 
reduced to the following form: 



Pp 



(Ax 2 + A 2 2 + A 3 2 ) V 



(16) 



other hand, as can be seen from (15), pp can ex- 
hibit non-trivial behaviors if the fluid is allowed to 
be anisotropic. In other words, in the presence of 
an anisotropic fluid, we may have cosmological models 
with non-trivial isotropization histories, e.g., models in 
which pp does not diverge as V — > and/or V — > +oo. 

2.4 Conservation laws for the anisotropic fluid and 
vacuum energy 

Allowing for anisotropy in the pressure of the fluid, and 
thus in its EoS parameter, gives rise to new possibilities 
for the evolution of the energy source. To see this, we 
first parametrize the energy-momentum tensor given in 
(2) as follows: 



= diag[p, -p x , ~p y , -p z ] 
= diag[l, -w x , -w v , -w z ]p 
diag[l, -w, -(w + 7), -(w + 6)]p, 



(17) 



where p x , p y and p z are the pressures and w x , w y and 
w z are the directional EoS parameters on the x, y and 
z axes, respectively; w is the deviation-free EoS pa- 
rameter of the fluid. The deviation from isotropy is 
parametrized by setting w x = w and then introducing 
skewness parameters S and 7, which are the deviations 
from w, respectively, on the y and z axes, w, 8 and 7 
are not necessarily constants and might be functions of 
the cosmic time, t. 

The conservation of the energy- momentum tensor of 
the fluid parametrized in (17), i.e., T^ v -, v — 0, leads to 
the following equation: 



p + 3(1 + w)pH + SpH y + ~fpH z = 0, 



(18) 



which is equivalent to the one defined by Barrow fc Turner^ 
( 19811 ) in case of an isotropic fluid. According to this 
equation, (16), pp decreases monotonically as the vol- 
ume of the universe increases and converges to null 
as V — > +00, while it diverges as V — > 0. On the 3(1 + w)H + SH y + jH z = 0, 



where the last two terms, the terms with 5 and 7, arise 
due to the anisotropy of the fluid. Now we can briefly 
examine what possibilities arise due to these two terms 
when the conventional vacuum energy (which is math- 
ematically equivalent to the cosmological constant A 
and which can be represented with an EoS parameter 
in the form of p = —p) is minimally altered in a way to 
wield an anisotropic EoS parameter. Obviously, when 
we consider the conventional vacuum energy, which is 
isotropic (i.e., 7 = £ = 0), if w — — 1, then the energy 
density (p) is necessarily constant and vice versa. How- 
ever, this is not the case when the energy-momentum 
tensor of the fluid is generalized to the form given in 
(17), namely, when an anisotropic EoS parameter is al- 
lowed for the fluid. In that case, if the energy density 
is assumed to be constant, p = const., (18) is reduced 
to the following equation: 



(19) 
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thus, w is not necessarily constant. Similarly, if w — 
— 1, (18) is reduced to the following equation: 



p + SpH y + 1P H Z = 0, 



(20) 



thus, the energy density is not necessarily constant ei- 
ther. 

2.5 Generalized Friedmann equation and de Sitter 
expansion 

Using the energy-momentum tensor parametrized in 
(17), the anisotropy energy density can be written in 
terms of the directional EoS parameters as follows: 



PP 



1 3 



Xj + / (wj — WijpVdt 



V and i > j. 



(21) 



Note that the behavior of the anisotropy energy den- 
sity is not only determined by the Xj constants and the 
volume of the universe as in (16), but also by the direc- 
tional EoS parameters (w x , w y and w z ) and the energy 
density of the fluid (p) via the integral term. 

We may write down the generalized Friedmann equa- 
tion by considering (13) and (14) and define the effec- 
tive energy density p e f , which determines the volumetric 
expansion rate of the universe, 



3H = p + pp = p ef . 



(22) 



From the definition of the mean Hubble parameter 
given in (9) the volume of the universe can be obtained 
as follows: 



V = Cl e^ Hd \ 



(23) 



where c\ is a positive constant of integration. Thus, 
using (22) in (23) the volume of the universe can be 
obtained in terms of the effective energy density, 



V = cie 



v3 / y/p c {dt 



(24) 



According to this, the condition for the de Sitter vol- 
umetric expansion does not correspond to a constant 
energy density of the fluid but to a constant effective 
energy density, i.e., p e f = P + Pp = const., which leads 
to 



V = cie 



/3p cf f 



(25) 



Obviously, the universe exhibits de Sitter expansion 
when pp is null in the presence of a positive cosmo- 
logical constant, i.e., in the presence of conventional 
vacuum energy that wields a constant energy density. 



However, if pp is not null, then p e f is not constant, but 
a decreasing function of the cosmic t as long as the uni- 
verse exp a nds, s ince pp oc V~ 2 . Thus, as also shown by 
iBeeshaml (jl994t ). the behavior of Bianchi type-I infla- 
tionary solutions cannot be of the pure de Sitter type, 
even when p = —p, but are of power-law type. 

What if we consider the anisotropy energy density 
together with a mixture of the conventional vacuum 
energy and perfect fluids? The conventional perfect flu- 
ids, i.e., radiation, pressureless matter etc. can be de- 
scribed by an EoS parameter in the form of p = wp and 
their energy densities change as V~^ 1+W K Thus, the 
cosmological constant will be dominant as V — > +oo, 
provided that w > — 1, and according to the cosm ic no- 
hair t heorem for Einstein gravity introduced by IWaldl 
(1983). the universe will exponentially evolve toward 
the de Sitter universe. In other words, Bianchi type-I 
models in the presence of a positive cosmological con- 
stant isotropize and their volumetric expansion rates 
approach de Sitter expansion as V — > +oo. On the 
other hand, according to (16), no matter how small the 
anisotropy energy density is, compared with the other 
sources, pp will eventually dominate any perfect fluid 
and govern the dynamics of the expansion in the very 
early evolution of the universe as V —> 0, provided that 
w < 1, i.e., the universe will approximate the Kasner 
vacuum solution. However, these cases are not implied 
once the implicitly assumed isotropy of the vacuum en- 
ergy is relaxed. This is because, according to (21), pp 
does not have to be proportional with V~ 2 in case of 
an anisotropic fluid, and it may not diverge as V — > 
and/or V — > +oo. 

Our final remark in this section concerns the con- 
dition for isotropization of a Bianchi type-I universe 
that exhibits de Sitter volumetric expansion. The con- 
dition for isotropization mentioned in Sect. 2.3 can be 
rewritten as follows: pp — > as t — > +oo since H is 
constant, and this leads to pp < as t — > +oo, since 
pp > by definition. This, in addition, leads to p > 
as t — > +oo from pp + p = 0. Thus, the condition for 
isotropization in the Bianchi type-I models that exhibit 
de Sitter volumetric expansion imposes on the fluid the 
constraint to behave like a phantom energy, i.e., to ex- 
hibit an increasing energy density as V increases. On 
the other hand, if the energy density of the fluid is con- 
stant (p = const.), then the anisotropy energy density 
is also constant (pp = const.). If p < 0, then pp > 0; 
thus for models in which the energy density of the fluid 
is decreasing, the anisotropy energy density (thus, the 
anisotropy of the expansion) increases as t increases. 
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3 Exponentially expanding LRS models in the 
presence of anisotropic fluid 

In the following, we examine the cosmological models 
that exhibit de Sitter volumetric expansion within the 
LRS Bianchi type-I framework and present two exact 
models. 

The LRS Bianchi type-I metric reads 



ds 2 = dt 2 - A{t) 2 dx 2 - B(t) 2 (dy 2 + dz 2 ). 



(26) 



Thus, H y = H z , and in the following they are repre- 
sented by H VtZ . The energy- momentum tensor given in 
(17) can be customized for the LRS Bianchi-I metric by 
choosing £ = 7, 



= diag[l, -to, -(w + 7), -(w + i)]p. 



(27) 



Considering (26) and (27), the Einstein field equations 
(4)-(7) can be reduced to the following system of equa- 
tions: 



B 2 AB 

Bi- + 2 AB =P > 



B 2 



B 
'B 



-wp, 



B BA A 

B + BA + A=- {W + 1)P - 



(28) 



(29) 



(30) 



Then we have initially five variables (A, B, p, w, 7) 
and three linearly independent equations, namely the 
three Einstein field equations (28)- (30). Thus we will 
need two additional constraints to close the system of 
equations. However, before introducing the constraints 
to close the system of equations, we can give p, w and 7 
in terms of the mean Hubble parameter and the direc- 
tional Hubble parameter on the x axis, manipulating 
the field equations (28)- (30), 



P = 3H 2 - - (H x - HY 



w 



A(H X ~3H)-1(H X -3H) 2 
m 2 -\ (H x - H) 2 



7 = 2jiL 



3 f t (H-H x ) + UH-H x )H 



(31) 



(32) 



(33) 



'AH 2 -l(H x - H) 2 

The anisotropy energy density can also be given in 
terms of H and H x by using (13), (14) and (28), 



PP = - A {H X -H) 2 . 



(34) 



One may check that the summation of (31) and (34) 
leads to p e { = p + pp = 3H 2 . 

As the first constraint to close the system of equa- 
tions, the effective energy density is assumed to be con- 
stant throughout the history of the universe: 

Pef = 3fc 2 , (35) 

where k is a positive constant and thus from (22) 

H = k, (36) 

which corresponds to the well-known de Sitter volumet- 
ric expansion, i.e., 



V = AB 2 = cie 



3fet 



(37) 



The solution of (33) by considering (31) and (36) gives 
the following equation for the evolution of the direc- 
tional scale factor on the x axis: 

A = Ke kt +& e ~ 3kt +&( e ~ 3kt f e 3kt r(t)dt-jr(t)dt) ^ ( 38 ) 

where k > and A are real constants and T(t) = jp 
is the skewness of the pressure. Using (38) in (37) the 
directional scale factor on the y and z axes is obtained 
as follows: 



B 



= (?iy /2 e kt-&e- 3kt -£(e- 3kt fe 3kt r(t)dt-fr(t) 



(39) 



Using the scale factors (38) and (39), the directional 
Hubble parameters arc obtained as follows: 



H x = k — Xe- 3kt - ^e- 3kt J r(t)e 3kt dt, 



k + -e 



-3kt 



+ 



e 3kt dt. 



(40) 



(41) 



Using (36), (40) and (41) in (10) the anisotropy of the 
expansion, A, and thus the anisotropy energy density, 
pp, are obtained as 

A _2p,_ 1 e - 6fct (3A + 2/ e ^r(t)^) 2 ^ 



3fc 2 18 



k 2 



It can be seen that the skewness of the pressure L(t) 
also contributes to the evolution of the cosmological 
parameters. 

One can immediately obtain the cosmological pa- 
rameters once a function for T(t) is chosen. Alterna- 
tively, one can obtain the parameters by introducing 
one more constraint on one of the cosmological param- 
eters; for instance, on the energy density of the fluid p 
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or on the EoS parameter. Trivially, one may also as- 
sume that the fluid is a perfect fluid (i.e., 7 = 0; thus 
r is null) and would th en obtain the solutions given by 
Kumar fc Singhl (|2007l) for an exponentially expanding 
Bianchi type-I universe. One can observe that for any 
perfect fluid, i.e., in case of T = 0, the initial anisotropy 
of the expansion dies away monotonically as t increases, 
while it may exhibit various non-trivial behaviors ac- 
cording to the function chosen for Tit), 

3.1 Model for p c { = const, and w = —I 

Note that we need one more constraint in addition to 
the one given in (35) to fully determine the system. 
Since the first assumption given by (35) causes de Sitter 
volumetric expansion, it seems reasonable to choose our 
last assumption as 



W 



-1. 



(43) 



If we solve the system of equations which consists of 
the three Einstein field equations (28)-(30) and the two 
constraints given by (35) and (43), we obtain the fol- 
lowing exact expressions for the scale factors: 



A 



Cl e kt Uk^+SXe- 3 ^^ 3 



B = e kt ink' 1 + 3\e 



3kt 



1/3 



(44) 



(45) 



The directional Hubble parameters are obtained as fol- 
lows: 



H>r — k 



6Afc 2 



Ke 3kt _|_ 3 \£ ' 



(46) 



H lt — k 



3Xk 2 



Ke 3kt _j_ 3 \£ 

The anisotropy parameter of the expansion and thus the 
anisotropy energy density are found to be dynamical, 



18X 2 k 



lPl = 

3 k 2 (ree 3fc * 



2 7,2 



3Afc)' 



(48) 



The energy density of the fluid is also found to be dy- 
namical in a way so as to secure p + pp = 3/c 2 , 



One can check that this solution yields (20). The 
positiveness condition on the energy density of the fluid 
(p > 0) imposes the restriction that the anisotropy of 
the expansion is smaller than 2 (A < 2) and that im- 
poses A to be positive. Thus, only those models with 
positive A values are viable, and hence in the following 
we will consider only the positive A values. 

A decreases monotonically as t increases and A — > 
as t -> +oo. Thus, the space approaches isotropy as 
t — > +oo in this model. On the other hand, A — » 
2 [I - k/(k + 3Afc)] 2 as t -> 0. 

p — » 3k 2 and 7 — > as t — > +00. That is, the EoS pa- 
rameter of the fluid isotropizes as t increases and mim- 
ics the cosmological constant ast^ +00. On the other 



27Affcj_ 



hand, p -> 3k 2 (l - (l + ^) 2 ) and 7 
which is always a negative number, as t — > 0. Therefore, 
as expected from the isotropization condition discussed 
in Sect. 2.5., p behaves like a phantom energy and in- 
creases as V increases. This result is consistent with 
the following situation: 7 is always negative and thus 
EoS parameters on the y and z axes are passing the 
phantom divide line, i.e., w y . z < — 1. 

3.2 Model for p c f — const, and p = const. 

While w is assumed to be —1 so as to close the system 
in the preceding section, it is allowed to be a function 
of the cosmic time t in this section, but, this time, the 
energy density of the fluid is assumed to be constant: 



p(t) = const. = p 



(51) 



If we solve the system of equations which consists of 
the three Einstein field equations (28)- (30) and the two 
constraints given by (35) and (51), we obtain the fol- 
lowing exact expressions for the scale factors: 



(47) A = K- 2 cie fet± tV9* 2 -3p* ) 



_ kt T ^9k 2 ^3pt 



B = ne 



(52) 



(53) 



The directional Hubble parameters are obtained as fol- 
lows: 



H x = k±-^/9k 2 -3p, 



(54) 



p = 3fc 2 1 



9X 2 k 



■2,1.1 



(ne 3kt + 3Afc)" 



(49) 



and the skewness parameter of the EoS parameter is 
obtained as follows: 



27A 2 fc 2 



Ke 3kt ( Ke 3kt _|_ 



(50) 



(55) 



The anisotropy parameter of the expansion and thus 
the anisotropy energy density are found to be non- 
dynamical, 



A = -Bl =2- 

3 k 2 3 k 2 ' 



(56) 
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The deviation-free EoS parameter of the fluid is ob- 
tained as follows: 

w= p-6g ±2 fc y ^-^ = _ i _A T ^ (g7) 

P P 1 2 

The skewness parameter of the EoS parameter is ob- 
tained as follows: 

7 = T3-v / 9fc 2 -3p = T-^. (58) 

^ 1 2 

One can check that this solution yields (19). The di- 
rectional Hubble parameters, the anisotropy of the ex- 
pansion (thus, the anisotropy energy density) and the 
directional EoS parameters are all constants through- 
out the history of the universe. 

It is obvious that p < 3k 2 , thus pp < 3k 2 and 
A < 2. We recover the conventional vacuum energy and 
isotropic expansion when p = 3k 2 , because w = — 1, 
7 = and A = in that case. However, when p < 3k 2 , 
both the expansion and the fluid deviate from isotropy, 
because w > —1, 7 < or w < — 1, 7 > and A > 
in this case. It can be observed that, while the EoS 
parameter on the x axis is in the quintessence region 
(w > —1), the ones on the y and z axes are in the 
phantom region (if + 7 < —1), or vice versa. According 
to this, while the expansion of the x axis acts so as to 
decrease the energy density of the fluid, the expansion 
of the yz plane acts so as to increase the energy density 
of the fluid or vice versa. However, in total, the decre- 
ments and increments compensate for each other and 
the energy density of the fluid does not change. 

4 Conclusion 

On motivating from the increasing evidence for the need 
of a geometry that resembles Bianchi morphology to 
explain the observed anisotropy in the WMAP data, 
we have discussed some features of the Bianchi type- 
I universes in the presence of a fluid that wields an 
anisotropic equation of state (EoS) parameter in gen- 
eral relativity. We have focused on those models that 
exhibit de Sitter volumetric expansion in the presence 
of a hypothetical fluid obtained by distorting the EoS 
parameter of the conventional vacuum energy in a way 
so as to wield anisotropy. 

We have also given two exact solutions within the lo- 
cally rotationally symmetric Bianchi type-I framework. 
In both models the effective energy density (the sum 
of the energy density of the fluid and the anisotropy 
energy density) has been assumed to be constant so as 
to secure the de Sitter volumetric expansion. 



In the first model, the directional EoS parameter on 
the x axis has been assumed to be -1. The anisotropy of 
the expansion, the energy density and the anisotropy of 
the fluid have been found to be dynamical. While the 
anisotropy of the expansion and the anisotropy of the 
fluid decrease and tend to null as the universe expands, 
the energy density of the fluid increases and approaches 
to its maximum value. The fluid approximates the con- 
ventional vacuum energy as the universe evolves. 

In the second model, the energy density of the fluid 
has been assumed to be constant. The anisotropy of 
the expansion and the anisotropy of the fluid have been 
found to be non-dynamical. When the maximum value 
of the energy density of the fluid (3k 2 , where A; is a 
positive constant) is considered, the universe expands 
isotropically and the fluid mimics the conventional vac- 
uum energy. Lower values of the energy density of the 
fluid give rise to an anisotropy both in the expansion 
and EoS parameter of the fluid. 
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